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Abstract. Our work focuses on the extended orbit problem: determin- 
ing whether a target vector space V may be reached from a point x under 
repeated application of a linear transformation A. We give a polynomial- 
time algorithm for the case when V is one-dimensional. We also show 
that when V is two- or three-dimensional, the problem is in the com- 
plexity class NP EqSLP , that is, NP with an oracle for deciding zeroness of 
arithmetic circuits, and hence in NP RP . 

1 Introduction 

The orbit problem is defined as follows: 

Given a rational matrix A and rational vectors x, y, decide whether there 
exists a non-negative integer n such that A n x = y. 

It may be interpreted as a reachability problem from point to point in a linear 
system with transitions described by A. In a well-known paper [T|, Kannan and 
Lipton showed a polynomial-time algorithm for the problem. The algorithm first 
reduces the orbit problem to the matrix power problem: 

Given a rational matrix A and a polynomial p £ Q[x], decide whether 
there exists a non-negative integer n such that A n — p(A). 

Then this problem is solved using the correspondence 

A n = p (A) a n = p (a) for all eigenvalues a of A 

Kannan and Lipton proved that if A has some eigenvalue which is not a root 
of unity, a polynomial bound m exists such that if A n = p(A), then n < m. 
With some effort they extended the decision method to cover matrices whose 
eigenvalues are all roots of unity, also in polynomial time. 

In their conclusion, they considered the following extended orbit problem: 

Given a rational matrix A and a vector subspace V specified by a basis 
of rational vectors, decide whether there exists a non-negative integer n 
such that A n x £ V. 

The extended orbit problem replaces the point y with a target space V. Kannan 
and Lipton conjectured that a polynomial-time algorithm exists for the case 
where dim(V) = 1, and suggested that decidability be investigated in the cases 



dim (V) = 2 and dim (V) = 3. Later work by Arvind and Vijayaraghavan [2] 
placed the original orbit problem in the logspace counting hierarchy GapLH. To 
the best of our knowledge, there has been no progress on the decidability of the 
extended orbit problem or its fixed-dimension cases. 

A related problem is the inhomogeneous extended orbit problem, where the 
target is an affine space V + z. The case dim (V) — is Kannan and Lipton's 
original problem. In general, the fc-dimensional inhomogeneous version reduces 
to the (k + l)-dimensional homogeneous version with the addition of a dummy 
component to all vectors. 

Another related problem is the chamber hitting problem. It replaces the target 
space with an intersection of half-spaces. In reference [3], the chamber hitting 
problem is related to the Pu-realizability problem on regular languages: deter- 
mine whether the language of a given finite deterministic automaton contains a 
word of the form #wi#W2# ■ ■ ■ #iJJN#, where N = 2™ is a fixed power of 2 and 
the words w%, . . . , wn are all the binary words of length n in some order. 

The orbit problem is also closely linked to Skolem's problem. Given a linear 
recurrent sequence described by its initial conditions and a recurrence formula, 
the problem is to decide whether is an element of the sequence. An equivalent 
formulation is to determine whether there exists a non-negative integer n such 
that x T A n y = for given vectors x, y and a square matrix A. The equivalence of 
these formulations is easy to show using the companion matrix of the sequence 
and the Cayley-Hamilton theorem. The dimension of A in the second formulation 
equals to the depth of the recurrence in the first. 

The decidability of Skolem's problem is open. A well-known result, the Skolem- 
Mahler-Lech theorem, states that the set of zeroes of any linear recurrence is the 
union of a finite set and finitely many arithmetic progressions |4I5I6I7| . More- 
over, it is known how to effectively compute the arithmetic progressions [7]. The 
main difficulty to deciding Skolem's problem is to determine whether the finite 
component of the set of witnesses is empty. 

Skolem's problem is trivial for recurrences of depth 1. A proof of decidability 
for depth 2 may be found in [8 . The problem was open for depth 3 and 4 until 
Vereshchagin [Q] and Mignotte, Shorey and Tijdeman [TO] proved it decidable in 
1985 using Baker's results for linear forms in logarithms [UJ. The case of depth 5 
was considered in Skolem's problem for arbitrary depth was shown NP-hard 
by Blondel and Portier [T3"] . 

Kannan and Lipton [JJ pointed out that Skolem's problem is easy to reduce 
to the orbit problem. Deciding whether x T A n y — for some n is equivalent to 
deciding whether A n y is in (span{x}) ± for some n. A reduction also exists in 
the other direction. Let (A, x, V) be an instance of the orbit problem, where V 
has basis {vx, . . . , Vf.}- Let {vk+i, ■ ■ ■ , v m } be a basis for V^. Then we have 



Given a yes-no oracle for Skolem's problem, the set of witnesses of a Skolem in- 
stance may be effectively computed. This is done by extracting the periodic sets 
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of zeroes and then iterating the sequence, using the oracle at each step to deter- 
mine if all the zeroes making up the finite component have been found. So we 
calculate the sets of witnesses for each instance and check whether the intersec- 
tion is non-empty. Or alternatively, if Sj are the linear recurrences corresponding 
to the Skolem instances vfA n x = 0, we can define the sequence 

m 

s(j)= E (*tf)) 2 
j=fe+i 

It is classical that linear recurrences are closed under addition and termwise 
multiplication, so S is a linear recurrence whose zeroes are exactly the witnesses 
to the orbit problem instance. 

Observe that in this reduction, the resulting Skolem instances have depth 
dependent on the size of the matrix in the orbit instance, not on the dimension 
of the target space. Since we are unable to solve Skolem instances for depth 
greater than 4, this severely limits the usefulness of the reduction. 

In this paper, we how to solve the orbit problem with a matrix of arbitrary 
size and target space of dimension k via Skolem's problem for depth k + 1. 
In particular, we show that the one-dimensional orbit problem is decidable in 
polynomial time and that the two- and three-dimensional versions are in the 
complexity class NP EqSLP and hence in NP RP . We first use a reduction to higher- 
dimensional versions of the matrix power problem, in the style of Kannan and 
Lipton. Then we construct an equivalent system of equations from the eigenval- 
ues of the given matrix. From the system, we extract tuples of equations and 
then apply known techniques from Skolem's problem to bound the exponent n. 

The structure of the paper is as follows. In section 2 we reduce the fixed- 
dimension orbit problem to the fixed-dimension matrix power problem. In section 
3, we prove that the one-dimensional version is decidable in polynomial time. In 
sections 4 and 5 we show decidability for the two- and three-dimensional cases. 
Finally, section 6 concludes the paper. Appendix A briefly outlines how algebraic 
numbers may be represented and used in arithmetic. Appendix B states some 
basic results from algebraic number theory. Appendix C states Baker's theorem, 
van der Poorten's theorem and some applications relevant to both the orbit 
problem and Skolem's problem. Appendices D, E and F give all the technical 
lemmas necessary for proving decidability of Skolem's problem in depth up to 4. 

2 Reduction 

2.1 Matrix power problem 

Suppose we are given a rational matrix A, a rational vector x and a target 
vector space V specified by a basis of rational vectors y 1 , . . . 7 y k . We wish to 
decide whether there exists n £ N such that A n x e V. 

Observe we can rescale A in polynomial time by the least common multiple of 
all denominators appearing in A. This reduces the problem to a version of itself 



but with the additional assumption that A is an integer matrix. This ensures 
the eigenvalues of A are all algebraic integers. 
Calculate 

v = max{m | x, Ax, . . . , A m x are linearly independent} 

Let B — {x, Ax, . . . , A v x}, U — span (B) and D — [x Ax . . . A v x\. It is clear 
that U is invariant under the linear transformation A, so consider the restriction 
of A to U. If [b , . . . , b u } T are the coordinates of A v+1 x with respect to B, then 
this restriction is described by the matrix 
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That is, if z = Bz B , then Az = DMzg. By induction, for all n e N, A n x = 
DM n XB, where x® = [1, 0, . . . , 0] T . Next we calculate a basis for W = U n V, 
let this basis be {w\, . . . , w t } and let Wi = Dw\ for all i. Now, 

A n x <E V A n x e W <=> M n x B e span {w[, w' t } 

Thus, the problem is reduced to a version of itself where the underlying matrix 
is invertible. Notice also that the matrix M describes a restriction of the linear 
transformation A, so its eigenvalues are a subset of the eigenvalues of A. In 
particular, since A was rescaled to an integer matrix, the eigenvalues of M are 
algebraic integers as well. 

Define the matrices Ti , . . . , T t by 

T i= [w[ Mw[...M v w[] 

We will show that M u xb G span {wi , . . . , w' t } if and only if M n g span {T 1 ,...,T t }. 
If for some coefficients aj we have 

t 

i=0 

then considering the first column of both sides, we have 

t 

M n x B = a * w 't 

Conversely, suppose M n xs — Y^\=a a i w 'i- Then note that xb, Mxb ■ ■ ■ , M v xs 
are just the unit vectors of size v + 1. Multiplying by M J for j = 0, . . . , v gives 
M n+ ^XB = J2l=o a i^ w \- The left-hand side is exactly the (j + l)-th column of 



M n , whereas M J u^ on the right-hand side is exactly the (j + l)-th column of 
T z . So we have M n = £* =0 a^i- 

Thus, we have reduced the orbit problem to the matrix power problem: 
determining whether some power of a given matrix lies inside a given vector 
space of matrices. Now we will perform a further reduction step. It is clear 
that out of the space T = span{T\, . . . ,T t } it suffices to consider only matri- 
ces of the shape p (M) where p is a polynomial. We find a basis for the space 
P = {p (M) \p£Q[x}} and then a basis { Pl (M) ,...,p s (M)} for PDT. Then 
M n £ T <^=> Al n £ P n T. We call this the polynomial version of the matrix 
power problem. Observe that dim (V) > dim(T) > dim(Tn P), so the dimen- 
sion of the target vector space does not grow during the described reductions. 
All described operations may be performed in polynomial time using standard 
techniques from linear algebra. 

2.2 Singular and nonsingular problem instances 

An instance (A, x, V) of the orbit problem is defined as non-singular if no quo- 
tient of two eigenvalues of A is a root of unity. In general, it is possible to reduce 
an arbitrary orbit problem instance to a set of non-singular instances as follows. 

Let L be the least common multiple of the orders of the quotients which are 
roots of unity. For each i £ {0, . . . ,L — 1}, we consider separately the problem 
of deciding whether there exists n £ N such that (A L ) (A l x) £ V. The original 
problem instance is positive if and only if at least one of these L instances is 
positive. 

It is easy to see that these instances are all non-singular. The eigenvalues of 
A L are exactly Af where A, are the eigenvalues of A. If for any two distinct such 
eigenvalues, say Xf ^ Xf, we have (Xf /Aj') = 1, then Xi/Xj must also be a root 
of unity. Then by the definition of L, Xf /Xf = 1, which gives the contradiction 



Observe that the reduction to the matrix power problem in the previous 
section preserves non-singularity, since the eigenvalues of the invertible matrix 
M are a subset of the eigenvalues of A. However, the non-singularity comes at 
the cost of higher complexity, as L may be exponentially large in the size of A. 
A deterministic machine would have to examine L problem instances to decide 
the original instance. Even if we allow non-determinism, we must still calculate 
A L , whose entries require space exponential in the size of the original problem 
instance. 

We show this reduction to make the point that if one only aims to prove 
decidability for the orbit problem, without regard for complexity, non-singularity 
of the problem instance may be assumed without loss of generality. Below we 
give tighter complexity upper bounds for the orbit problem with a target space of 
dimension at most 3 and explicitly handle singular problem instances, but if we 
only aimed at proving decidability, the argument could be shortened significantly. 

Finally, we point out that a similar trick was used by Vereshchagin [5] to 
prove decidability of Skolem's problem for recurrences of depth 3 and 4. Given 




a singular recurrent sequence u (n) , one may take 

L = 1cm {m | Xi/Xj is an m-th root of unity} 

where Xi are the roots of the characteristic polynomial of the sequence. Then 
each sequence Vi (n) = u (Ln + i) for i G {0, . . . , L — 1} is non-singular in the 
sense that dividing two of its roots does not give a root of unity. 

2.3 Towards a system of equations 

Suppose now we have an instance (A,p%, . . . ,p s ) of the polynomial version of 
the matrix power problem. Calculate the minimal polynomial Ja (x) of A and 
obtain canonical representations of its roots ai, . . . , ctk, that is, the eigenvalues 
of A. This may be done in polynomial time, see Appendix A. Throughout this 
paper, for an eigenvalue a, we will denote by mul (on) the multiplicity of on in 
the minimal polynomial of the matrix. 

Fix an exponent n € N and coefficients a\, . . . ,a s € C and define the poly- 
nomials P (x) — a iPi ( x ) an d Q (x) = x n . It is easy to see that 

Q{A) = P{A) 

if and only if 

v ? : € {1, . . . , k} .Vj e {0, . . . , mul (o<) - 1} .P {j) (a t ) = qCj) ( ai ) (1) 

Indeed, P — Q is zero at A if and only if Ja (x) divides P — Q, that is, each ai is 
a root of P — Q with multiplicity at least trail (a^). This is equivalent to saying 
that each is a root of P — Q and its first mul (ai) — 1 derivatives. 

Thus, in order to decide whether there exists an exponent n and coefficients 
ai such that A n = J2t=i a iVi {A), it is sufficient to solve a system of equations 
(JTJ where the unknowns are n and ai, . . . ,a s . Each eigenvalue ai contributes 
mul (ai) equations which specify that P (x) and its first mul (ai) — 1 derivatives 
all vanish at a^. 

For brevity in what follows, we will denote by eq (ai,j) the j'-th derivative 
equation contributed to the system by ai, that is, P^' (ai) = Q^) (ai). This 
notation is meaningful only for < j < mul (ai). We will also denote by Eq (ai) 
the set of equations contributed by ai to the system: 

Eq (a^ = {eq (oj,0) , . . . , eg (a,-, mul (on) - 1)} 

For example, if j a (x) has roots a\, 02, 03 with multiplicities mul (ai) — i and 
the target space is span {p± (A) ,p2 (A)} then the system contains six equations: 

a-iPi (at) + a 2 p 2 ("l) 
a-iPi (a 2 ) + a 2 p2 («a) 
aip[ (a 2 ) + a 2 p 2 
a-iPi ("3) + 02^2 ("3) 
aiPi ("3) + 02^2 ("3) 
aiPi ("3) + a2P 2 ' ( a 3) 



n (n — 1) d!i~ 



Then eq (03, 0) is the equation 



«3 = ( a 3) + a2P2 (a 3 ) 
and Eq (a 2 ) is the two equations 

ag = aipi (a 2 ) + a 2 P2 ("2) 

na^' 1 — a iPi ( a 2) + ( a z) 

3 One-dimensional version 

Suppose we are given a one-dimensional matrix power problem instance (A,p) 
and wish to decide whether A n 6 span {p (A)} for some n. We have constructed 
a system in the exponent n and the coefficient a as in ([lj. For example, if the 
roots of /a (x) are a±,a2, 013 with multiplicities mul (ccj) = i, the system is: 
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In this section we will describe how such systems may be solved in polynomial 
time. We allow the problem instance to be singular, that is, the ratios of eigen- 
values of A may be roots of unity. 

The strategy is to consider quotients of equations. This eliminates the un- 
known coefficient a and leaves only the exponent n. First, we perform some 
preliminary calculations. 

1. We check whether a = has a corresponding n which solves the matrix equa- 
tion A n — ap(A). This may be done using Kannan and Lipton's algorithm 
for the original orbit problem. If this is the case, we are done. Otherwise, 
assume o^O. 

2. Let c = maxj {mul (aj)}. We check for all n < c whether A n is a multiple of 
p (A). If so, we are done. Otherwise, assume n > c. 

3. We check whether on — for some i. If so, then all of the equations Eq (a;) 
are of the form = ap^ (0), which is equivalent to = p^ (0). We can easily 
check whether these equations are satisfied. If so, we dismiss them from the 
system without changing the set of solutions. If not, then we are done. Now 
we assume a^ ^ for all i. 

4. Finally, we check whether the right hand side apW (a,;) of some equation is 
equal to 0, using a polynomial division of p^ (x) by the minimal polynomial 
of OLi. If this is the case, then the problem instance is negative, because the 
left-hand sides are all non-zero. 



Let eq(ai,k) /eq(otj,t) denote the equation obtained by dividing eq(oti,k) by 
eq (ctj, t), that is, 



n (n - 1) . . . (n - k + 1) a"~ fe _ pW (a*) 
n (n - 1) . . . (n - t + 1) ~~ (^) 

We compute canonical representations of all quotients ctijotj. and consider three 
cases. 

Case 1. Some quotient ctijctj is not a root of unity. Then eq(pn,G) and 
eg (aj, 0) together imply eg (cti, 0) /eg (<x,-, 0), that is, 



<x,/ p(<x,-) 

By Lemma ^ in Appendix A, we can compute canonical representations of 
p(cti) /p(ctj) and cti/ctj in polynomial time. Then by Lemma[5]in Appendix C, 
n is bounded by a polynomial in the input. We check A n £ span {p (A)} for all 
n up to the bound and we are done. 

Case 2. All quotients otifoij are roots of unity, and all roots of Ja (x) are 
simple. Then the system is equivalent to 

a " a A e 9( a *>°) 
a = — -, — r A A — - 

Pip-x) l^eqyotj^) 
It is sufficient to determine whether there exists some n which satisfies 

A ^ (q»>0) (9] 
(\eq(a v 0) U 

Consider each equation eq (a^, 0) j eq (a^, 0): 

' ai\ n p{oti) 



OtjJ p(otj) 



(3) 



Suppose aj /ckj is an r-th root of unity. If the right-hand side of (|3]) is also an r-th 
root of unity, then the solutions of ([3| are n = t mod r for some t. If not, then 
Q has no solution, so the entire system ([TJ has no solution, and the problem 
instance is negative. By Lemma[l] we can determine in polynomial time whether 
the right-hand side of ^ is a root of unity, and if so, calculate t. We transform 
each equation in Q into an equivalent congruence in n. This gives a system 
of congruences in n which is equivalent to We solve it using the Chinese 
Remainder Theorem. The problem instance is positive if and only if the system 
of congruences has a solution. 

Case 3. All quotients oti/uj are roots of unity, and Ja (x) has repeated 
roots. We transform the system into an equivalent one in the following way. 
First, we include in the new system all the quotients of equations eq(ai,0) 



as in Case 2. Second, for each repeated root ai of /a(x), we take the quo- 
tients Aj"To^ Ql ' ) 2 e Q j) I e 1 { a i> 3 : + !)■ Third, we include the equation a = 
ai/p(ai). 

mul(oti)-2 ,. 

A £? a "° AA A qM a 8 = ^ 

We solve the first part of the system as in Case 2. If there is no solution, 
then we are done. Otherwise, the solution is some congruence n = t\ mod t 2 - 
For the remainder of the system, each ratio eq (ai, j) / eq (ai,j + 1) contributed 
by a repeated root ai has the shape 

(i) 



ai p 1 - 1 ' (ai) 



j pO'+i) (a^ 



which is equivalent to 



P {j+1) (oti) 

For each such equation ([!]), we calculate the right-hand side and check whether 
it is in N. If not, then the system has no solution. Otherwise, Q points to a 
single candidate n . We do this for all equations where n appears outside the 
exponent. If they point to the same value of n, then the system is equivalent to 

n = t\ mod ti 

n = rio 

a = oty/q(ax) 

We check whether uq satisfies the congruence and we are done. 



4 Two-dimensional version 

Suppose we are given (A,pi,p 2 ) where A is a square matrix and Pi,P2 are 
polynomials. We wish to decide whether there exists n £ N such that A" E 
span {p\(A),p2(A)}. In this section we will show that this problem is in the 
complexity class NP EqSLP , and hence in NP RP . The instance is allowed to be sin- 
gular, that is, there may exist distinct eigenvalues of A whose ratio is a root of 
unity. We have derived a system of equations ([T]). 

The broad strategy will be to choose a tuple of equations from the system 
(JlJ, obtain a Skolem instance of depth 3 from this tuple and hence obtain a 
bound m on the exponent n such that if 

A n G span{p 1 (A) lP2 (A)} 

then n < m. This bound will be at most exponential in the size of the input, 
so an NP machine will be able to guess an exponent n up to the bound. Then 



the machine calculates A", representing numbers as arithmetic circuits, and 
expresses membership in the target vector space as an instance of the EqSLP 
problem: determining whether a given arithmetic circuit evaluates to 0. 

For example, suppose A has three distinct eigenvalues a,/?, 7, and consider 
the tuple of equations eq (a, 0), eq (/?, 0), eq (7, 0): 

P n = <*i [pi 08) + a 2 P2 {[3) 
T) \Pih)J \P2h)J 

If the vectors on the right-hand side are linearly independent, then this triple 
states that the point in A 3 described by the left-hand side lies on the plane in 
A 3 described by the right-hand side. We can calculate the normal (Ai, A2, Az) T 
of the plane to obtain 

A ia n + A 2 p n + A37" -0 

Now observe that the left-hand side as a function of n satisfies a linear recurrence 
of depth 3 over A. Provided that the ratios of a, /3, 7 are not roots of unity, 
Lemmas [6j [7] [8] give a bound on n which is at most exponential in the input 
size, as desired. Problems arise, however, when this linear recurrence is allowed 
to be singular. For example, suppose that A3 = 0, and let a/j3 and —A2/A1 be 
roots of unity of the same order. Then the zeroes of the recurrence are a full 
arithmetic progression, and this Skolem instance fails to give a bound on n. 

Similarly, if the matrix A has only two eigenvalues a, /3, but one of them, 
say a, is repeated, then we may consider the triple eq (a, 0), eq (a, 1), eq (/3, 0). 
Using similar reasoning, we see that 

A ia n + Azna 71 - 1 + A 3 (3 n = 

must hold for some effective algebraic constants Ai , A2 , A3. This corresponds to 
a Skolem instance of depth 3 where one of the characteristic roots is repeated. 
Now provided that a//3 is not a root of unity, we have an exponential bound on 
n from Lemma[9j However, if a//3 and — A3/A1 are roots of unity, and A2 = 0, 
then the Skolem instance could have infinitely many solutions, failing to give a 
bound on n. 

In order to explicitly handle the case of singular orbit instances, we will 
consider the relation ~ on the eigenvalues of A, defined by 

a ~ P if and only if a//3 is a root of unity 

It is clear that ~ is an equivalence relation. The equivalence classes Ci, . . . ,Ck of 
~ are of two kinds. First, a class can be its own image under complex conjugation: 

Ci — {a I a e Ci} 

Each such self-conjugate class . . . , a s } has the form {awi, . . . , au> s } where Wj 
are roots of unity, and \a.j\ = a € Rfl A. Call this a the stem of the equivalence 
class d. Second, if an equivalence class is not self-conjugate, then its image 



under complex conjugation must be another equivalence class of ~. Thus, the 
remaining equivalence classes of ~ are grouped into pairs (Cj , Cj ) such that Cj = 
{x | x € Cj} = Cj. In this case, we can write Cj and Cj as 

Ci = {Xloi, . . . , Xu> s } 

C 3 = {Acji, . . . , Aw s } 

where cjj are roots of unity and A € A with arg (A) /2-7T ^ Q. Call A the stem 
of Cj and A the stem of Cj. Observe that the stems of self-conjugate classes are 
distinct positive real numbers, and that the ratio A/A of the stems of paired 
classes cannot be a root of unity. Recall also that we can assume the eigenvalues 
of A are algebraic integers, as a by-product of the reduction from the orbit 
problem. Since roots of unity and their inverses are algebraic integers, it follows 
that the stems of equivalence classes must also be algebraic integers. 
Let 

Eg (C) = |J Eg (a) 

denote the set of equations contributed to the system by the eigenvalues in C, 
and let 

Eq(C,i)= (J {eg(a,i)} 

aeC 
i < mul [a] 

denote the set of i-th derivative equations contributed by the roots in C. 
We will case split on the equivalence classes of ~. 

Case I. Suppose ~ has exactly one equivalence class C = {awi,...,au s }, 
necessarily self-conjugate, with stem a. Let L be the least common multiple of 
the orders of oji, . . . ,w s . Since L is at most exponentially large in the size of the 
input, it can be written using polynomially bits, so our NP machine will guess 
r = n mod L. Now consider the set of equations Eg (C, 0): 

(au}{) = aip (qwi) -I- C12P (awj.) 

(auj s ) n = aip (aw s ) + a 2 p (auj s ) 

Having guessed r, we can easily calculate w™, . . . , cj" in polynomial time, since 
Mi are roots of unity whose order divides L. Then the equations Eq(C,0) are 
equivalent to 

(5) 

where B is an s x 2 matrix over A, computable in polynomial time. Next we 
subtract the first row of (|5| from rows 2, . . . , s, obtaining 

a n = tpxCLi + ip 2 a 2 A 
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Here ( ipi (f2) is the first row of the matrix B, and B' is the bottom s — 1 rows 
of B. Thus, Eq (C, 0) is equivalent to 



a n = (fiax + Lp 2 a 2 

together with the constraint that (ai a 2 ) T must lie in the nullspace of B' . We 
calculate the nullspace of B' directly. If its dimension is less than 2, then we 
have a linear constraint on ai,a 2 . This constraint is of the form = ka 2 when 
the nullspace of B' has dimension 1, and is a\ = a 2 = when the nullspace is 
of dimension 0. In both cases, we substitute into the system Q, and solve the 
resulting lower-dimensional system using the algorithms for the one-dimensional 
orbit problem and Kannan and Lipton's original orbit problem. In the case when 
the nullspace of B' has dimension 2, then the linear constraint is vacuous, and 
Eq (C, 0) is equivalent to a 11 = (piai + (p 2 a 2 . 

In the same way we collapse Eq (C, 1) into one first-derivative equation 

fia n_1 = (pzai + (p^a 2 

using our guess for n mod L. We do this for all Eq (C,i), obtaining a system 
of equations equivalent to based on the stem of C, rather than the actual 
eigenvalues in C. Denote the resulting set of equations by T (Eq(C)). 

If some eigenvalue x £ C has mul (x) > 3, then J- (Eq (C)) contains the 
following triple of equations: 



(0) 




n (n — 1) a 



If the vectors on the right-hand side of ^ are linearly independent, then they 
specify a plane in A 3 , and the triple states that the point on the left-hand side 
must lie on this plane. We calculate the normal (Ai A 2 A3 ) of the plane and 
obtain 

A x a n + A^ia 71 - 1 + A 3 n (n - 1) a"~ 2 = 

This is a quadratic equation in n. It has at most two solutions, both at most 
exponentially large in the size of the input, so we are done. If the vectors on the 
right-hand side of ^ are linearly dependent, then we may equivalently consider 

a n 

na n ~ 1 I = " 1 
n (n — 1) a n ~ 2 

We divide the first equation by the second to obtain 

°L — fl 
n ip 3 




which limits n at most one, exponentially large, candidate value. 



If all eigenvalues x in C have mul (x) < 2 and at least one has mul (x) = 2. 
then T (Eq (C)) consists of exactly two equations: 



If (tpi if3) T and (t/J2 V4) T are linearly independent, then the right-hand side 
of (I7| spans all of A 2 as ai, a 2 range over A. so the problem instance is trivially 
positive. Otherwise, we may equivalently consider 

a n \ _ ( Vl 

and limit n to at most one candidate value which is exponentially large in the 
input size. 

Finally, if all eigenvalues x in C have mul (x) = 1, then T (Eq (C)) contains 
only the equation 

a™ = anpi + a 2 ifi2 

which has a solution if and only if at least one of ipi , <f2 is non-zero. 

Case II. Suppose ~ has exactly two equivalence classes, C\ and C2, with 
respective stems a and f3, so that 

Ci = {out, . . . ,acj s } 

C 2 = {M,...,^} 

The classes could be self-conjugate, in which case a, (3 £ A n K, or they could 
be each other's image under complex conjugation, in which case a = j3. In both 
cases, a//3 is not a root of unity. 

As in the previous case, we define L to be the least common multiple of the 
orders of u>i, . . . ,uj s , u)[,. . . and our NP machine guesses r — n mod L. We 
transform the system Eq (C\) A Eq (C2) into the equivalent system T (Eq (Ci)) A 
F (Eq (C2)). If all eigenvalues x of A have mul (x) = 1, then the system consists 
of two equations, one for each equivalence class of ~: 

a n = ai^i + a 2 </52 
(3 n = ai^ 3 + a 2 <^4 

If ((fx V3) T and (tp2 l P4) T are linearly independent, then the problem instance 
is trivially positive. Otherwise, it suffices to look for n which satisfies 

a n = a±(pi 
P» = a^ 3 

and hence 

/ \ n 



A bound on n follows from Lemma [5] This argument relies crucially on the fact 
that a/ f3 is not a root of unity. 

If some eigenvalue x of A has mul(x) > 2, say x £ C\, then the system 
contains the following triple of equations: 



(8) 



If the vectors on the right-hand side of (|8| are linearly dependent, so that the 
right-hand side describes a space of dimension 1 , it suffices to look for solutions 
to 

ai I <^3 




Then dividing we obtain 



<^3 



which limits n to at most one, exponentially large candidate value. Otherwise, if 
the vectors on the right-hand side of ^ are linearly independent, we calculate 

the normal (A\ A 2 A 3 ^j T to the plane described by them and obtain 

A ia n + A 2 na n ~ l + A 3 p n = 

A bound on n which is exponential in the size of the input follows from Lemma 
[9] This relies on the fact that a/ (3 cannot be a root of unity. 

Case III. Suppose ~ has at least three equivalence classes. Then we can 
choose eigenvalues a, j3, 7, each from a distinct equivalence class. Then consider 
eq(a, 0), eq(/3,0) and eg (7,0): 

[Pl( a )\ (P2{0L)" 

= ai Pi 09) + a 2 \p 2 (P) 

\pi(i)J \P2(j), 

If the vectors on the right-hand side are linearly independent, we calculate the 
normal (Ai, A 2 , A^) T of the plane on the right-hand side to obtain 

A ia n + A 2 p n + A 3l n = 

This is a non-singular Skolem instance of depth 3, so a bound on n follows from 
Lemmas [6] [7] [8] If the vectors on the right-hand side are not linearly independent, 
then we may equivalently consider 

/ Pi 

= ai Ui 09) 

\Pi (7), 





which gives 

M"_ Pi (a) 
W Pi (/?) 

An exponential bound on n follows from Lemma [5] because a//3 is not a root of 
unity. 

Thus, we have now shown that in all cases, an NP machine may compute a 
bound m such that if 

A n x € span {y, z} 

then n < m. This bound is at most exponential in the size of the input. From 
here it is easy to argue membership in NP EqSLP . The machine guesses some n 
up to the bound, ensuring that this n is consistent with the guess for n mod L. 
Now we need to compute A n x and check whether it is in the target vector space. 
Since n is exponential in the size of the input, the entries of A n x are, in general, 
doubly-exponential. That is, they require an exponential number of bits to write 
down. However, the entries of A n x may easily be represented as polynomial- 
sized arithmetic circuits. We consider all projections of A n x, y and z to three 
coordinates, and for each projection we expression the question of linear inde- 
pendence as the zeroness of a 3 x 3 determinant, also expressed as an arithmetic 
circuit. It is clear that n is a witness to the problem instance if and only if for 
any projection to three coordinates, A n x, y and z are linearly dependent. This 
is easy to determine with an EqSLP oracle, so we have membership in NP EqSLP . 
It is known that EqSLP C coRP |14| . so we also have membership in NP RP . 

5 Three-dimensional version 

Suppose we have a problem instance (A,p\,p2,P3) and wish to decide whether 
A n s span {pi (A) ,p2 (A) ,p^ (A)} for some n. As before, we have constructed a 
system (JT|) in n and the coefficients a%, ct2, as- The eigenvalues of A are algebraic 
integers and do not include 0. In this section we will show that there exists an 
effective bound m which is at most exponentially large in the size of the input, 
such that 

A n G span { Pl {A) ,p 2 (A) ,p 3 (A)} =4> n < m 

Then by the same reasoning as in the two-dimensional case, we will have mem- 
bership in NP EqSLP and NP RP for the three-dimensional orbit problem. 

Following the strategy of the two-dimensional case, we will select tuples of 
equations and obtain a bound on n using the lemmas for Skolem's problem for 
recurrences of depth 4 in Appendix F. We will again case split on the equivalence 
classes of the relation ~. 

Case I. Suppose there are at least two pairs of classes (Cj,QJ, (Cj,Cj) which 
are not self-conjugate. Then let a G Cj, = a G Ci, 7 G Cj, S = 7 G Cj. Then we 



consider the tuple of equations 



( an \ 




(Pi (<*)\ 




( P2 (a) \ 




( P3 (a) \ 






Pi ifi) 


+ a 2 


P2 (P) 


+ 0*3 


P3 (P) 


7" 


= ai 


Pi (7) 


Pi (7) 


P3 (7) 


w 




\pi (S)J 




\P2 (S)J 




\P3 (6) J 



(9) 



If the vectors on the right-hand side are linearly dependent, then we rewrite the 
right-hand side as a linear combination of h < 3 vectors and obtain a bound on n 
as we did for tuples of equations in the one- and two-dimensional orbit problem. 
If the vectors on the right-hand side of (JoJ) are linearly independent, then we 
calculate the normal of their three-dimensional subspace of A 4 , obtaining an 
equation 

A x a n + A 2 p n + A 3l n + A 4 S n = (10) 

and hence an exponential bound on n from Lemmas [12] and |13| We are relying 
on the fact that the ratios of a, /3, 7, 8 are not roots of unity. Notice that we need 
(a, /3) and (7, 8) to be pairwise complex conjugates in order to apply Lemma 13 
Case II. Suppose now that there is exactly one pair of classes (Ci,Ci) which 
are not self-conjugate. In general, for any eigenvalue x of A we must have 
mul (x) = mul (x). Therefore, if any eigenvalue a e Cj has mul (a) > 1, we can 
select the non-singular tuple of equations eq (a, 0), eq (a, 1), eq (a, 0), eq (a, 1): 



/ a" \ 




(Pi (a)\ 




f P2 (") ^ 




f P3 (a) \ 


a™ 




Pi (a) 




P2 (a) 




P3 (a) 




= ai 


Pi (") 


+ a 2 


P2 ( a ) 


+ a 3 


P'-i (") 


{net 71 - 1 ) 




\P'i(a)J 




\P2(a)J 







This leads to a non-singular Skolem instance of depth 4 over A for a recurrence 
sequence with two repeated characteristic roots: 



Axet 71 + A 2 a n + A 3 na n 1 + A 4 na n ~ 







An exponential bound on n follows from Lemma 10 since a/ a is not a root of 
unity. 

We can now assume that eigenvalues in Ci and Ci contribute exactly one 
equation to the system. Now we perform the collapsing operation used in the 
two-dimensional orbit problem, transforming Eq (Cj) A Eq (C^j into the equiva- 
lent T (Eq (d)) A T (Eq (Cj)) . Since all eigenvalues in d and Ci contribute one 
equation each, T (Eq (Ci)) A J- (Eq (Ci)) is just 

A" = aupt + a 2 (fi2 + 03^3 
A = 01^4 + a 2 ip 5 + a 3 ip 6 

where A, A are the stems of Ci and Ci. We perform the collapsing operation to all 
self-conjugate classes as well, reducing the system of equations to an equivalent 
system based on the stems of the equivalence classes, not the actual eigenvalues 
of A. This is beneficial, because the stems cannot divide to give roots of unity, 



so we can use 4-tuples of equations to construct non-singular Skolem instances 
of order 4. 

If there are at least two self-conjugate equivalence classes, with respective 
stems a, (3, we take the tuple 

= a\ip\ + a 2 tp 2 + a-3f3 
= aiip4 + a 2 ip 5 + a 3 tp 6 
= Cli</>7 + a 2 ip$ + «3<y59 
ai^io + a 2 ipn + a 3 tpi 2 

and obtain the non-singular Skolem instance 

AiA n + A 2 T + A 3 a n + A A p n = 

Then we have a bound on n from Lemmas [12] and [13] Similarly, if there is only 
one self-conjugate equivalence class, with stem a, but some of its eigenvalues are 
repeated, we use the tuple 

ai<pi + a 2 ip 2 + a 3 tp 3 
aiifi + a 2 ip 5 + a 3 ip 6 
ai<P7 + a 2 ifs + a 3 (p 9 
aifw + a 2 ipn + a 3 ip X2 

to obtain the non-singular instance 

AiA" + A 2 X n + A 3 a n + A^na 11 - 1 = 

which gives a bound on n according to Lemma |11| If there is exactly one self- 
conjugate class, with stem a, containing no repeated roots, then the system 
consists of three equations: 

A" = aitpi + a 2 ip 2 + a 3 tp 3 
A = (nipi + a 2 (f5 + a>3¥>6 
a" = ai<^7 + a 2 ip$ + a 3 tp 9 

Depending on whether the vectors ( ipi ip^ tpi ) T , ( ip 2 ip$ ipg ) T , ( ip 3 ip§ ip 9 ) T 
are linearly independent, this is either a trivially positive instance, or a lower- 
dimensional non-singular instance. Finally, if there are no self-conjugate classes, 
the system consists of only two equations: 

A" = ax<fi + a 2 if 2 + a 3 Lp 3 

A = (1^4 + a 2 lf 5 + a 3 (p 6 

Again, depending on the dimension of 

this is either a trivially positive instance, or a lower-dimensional non-singular 
one. 



A" = 
A" = 
a n = 

/3 n = 



A" = 
A" = 
a n = 



Case III. All equivalence classes of ~ are self-conjugate. As above, we use the 
collapsing operation, and pick out 4-tuples of equations. We rely on the fact that 
stems of classes are distinct real algebraic numbers to ensure that the resulting 
Skolem instances are non-singular and give us the desired bound on n. 

6 Complexity and conclusion 

We have shown that the higher-dimensional orbit problem is decidable in poly- 
nomial time when dim (V) = 1. We have also shown membership in NP EqSLP in 
the cases dim (V) G {2, 3}. It is known [2] that EqSLP C coRP, so membership 
in NP RP follows immediately. 

The connection between Skolem's problem and the orbit problem in fixed di- 
mension is evident. A Skolem problem instance requires us to find the zeroes of a 
quasi-polynomial X)i=o a ?Pi whereas an orbit problem instance is equivalent 
to a system of equations from which we may extract quasi-polynomial equations. 

In general, given a Skolem instance 

3n.x T A n y = 

where A has size fc x fc, we can reduce to the equivalent orbit instance 

A n y G (span{x}) 1 ' 

which has a target vector space of dimension k — 1. In the other direction, 
starting with an orbit instance with a fc-dimensional target, we can construct an 
equivalent system of equations, extract (A; + l)-tuples, obtain a normal vector of 
a fc-dimensional subspace and hence an equation which corresponds to finding 
the zeroes of a recurrence of depth fc + 1. 

It is interesting to note that at first glance, the orbit problem in fixed dimen- 
sion appears more difficult that Skolem's problem for a fixed depth, due to its 
unbounded matrix. On the contrary, the presense of more eigenvalues allows us 
more freedom when choosing equations from which to extract a bound on the 
exponent n. Repetitions in the roots of Ja simplify matters greatly by bringing 
n into the base position. A Skolem problem instance offers no such freedom. 

Because of the similarities in the mathematics of these two problems, we 
believe they should be explored in conjunction. Given that Skolem's problem 
for depth 3 and 4 was only solved with the advent of Baker's theorem, it is 
unsurprising that we required the same results to prove decidability for the two- 
and three-dimensional orbit problem. 

A Algebraic numbers and operations on them 

A complex number a is algebraic if there exists a polynomial p G Q [x] such 
that p (a) = 0. The set of algebraic numbers, denoted by A, is a field under 
the usual arithmetic operations. The minimal polynomial of a is the unique 



monic polynomial of least degree which vanishes at a and is denoted by f a (x). 
The degree of a £ A is defined as the degree of its minimal polynomial and is 
denoted by n a . The height of a is defined as the maximum absolute value of a 
coefficient in the minimal polynomial of a and is denoted by H a . The roots of 
f a (x) (including a) are called the Galois conjugates of a. The absolute norm 
of a, denoted Af a bs (en), is the product of the Galois conjugates of a. By Viete's 
laws, we have 

Nabs (a) = (-1)"° - h 

where a, b are respectively the free term and the leading coefficient of f a (x). It 
follows that Nabs (a) € Q- An algebraic integer is an algebraic number a such 
that f a £ Z[x\. The set of algebraic integers, denoted Az, is a ring under the 
usual addition and multiplication. 

The canonical representation of an algebraic number a is its minimal poly- 
nomial fa (x), along with a numerical approximation of Re (a) and Im (a) of 
sufficient precision to distinguish a from its Galois conjugates. More precisely, 
we represent a by the tuple 

(f a ,x,y,R)e{Q[x]xQxQxQ) 

meaning that a is the unique root of f a inside the circle centred at (x, y) in the 
complex plane with radius R. A bound due to Mignotte [T5] states that for roots 
oii 7^ otj of a polynomial p (x), 

l a *- Q ji > n (n+l)/2ffn-l (H) 

where n and iJ are the degree and height of p, respectively. Thus, if R is restricted 
to be less than a quarter of the root separation bound, the representation is well- 
defined and allows for equality checking. Observe that given f a , the remaining 
data necessary to describe a is polynomial in the length of the input. It is known 
how to obtain polynomially many bits of the roots of any p £ Q [x] in polynomial 
time [IB] , 

From here onwards, when we say an algebraic number a is given, we will 
assume we have a canonical description of a. We will denote by ||a|| the length 
of this description, assuming that integers are expressed in binary and rationals 
are expressed as pairs of integers. Observe that \a\ is an exponentially large 
quantity in ||a|| whereas ln|a| is polynomially large. Notice also that l/ln|a| is 
at most exponentially large in ||a||. For a rational a, ||a|| is just the sum of the 
lengths of its numerator and denominator written in binary. For a polynomial 
p € Q [x], will denote X^ILn \\Pi\\ wnere n is the degree of the polynomial and 
Pi are its coefficients. 

Lemma 1. Given canonical representations of a, ft £ A and a polynomial p £ 
Q[x], it is possible to compute canonical descriptions of a ± /3, andp(a) 
in time polynomial in the length of the input (that is, in \\a\\ + \\/3\\ + \\p\\). 



Proof. The resultant of f a (x — y) and fp (y), interpreted as polynomials in y 
with coefficients in Q[x], is a polynomial in x which vanishes at a + (3. We 
compute it in polynomial time using the Sub-Resultant algorithm (see Algo- 
rithm 3.3.7 in |17| ) and factor it into irreducibles using the LLL algorithm [18 . 
Finally, we approximate the roots of each irreducible factor to identify the mini- 
mal polynomial of a + (3. The degree of a + /? is at mos t n a np, while its height is 



bounded by H a+ p < H™ a Hp |19j . Therefore, by (111, a polynomial number of 
bits suffices to describe a + f3 unambiguously. Similarly, we can compute canon- 
ical representations of a — f3, a(3 and a//3 in polynomial time using resultants, 
see HH]. 

To calculate p (a) we repeatedly use addition and multiplication. It suffices 
to prove that all intermediate results may be represented in polynomial space. It 
is clear that their degrees are at most n a , but it is not obvious how quickly the 
coefficients of their minimal polynomials grow. However, there is a simple reason 
why their representation is polynomially bounded. Let A be the companion 
matrix of f a . Then p (a) is an eigenvalue of p (A). We can calculate p (A) using 
only polynomial space. Then from the formula 

n 

det (XI -p(A))= £ S gn(o-)l[(\I-p(A)) lMl) 

it is evident that the coefficients of the characteristic polynomial of p (A) are 
exponentially large in the length of the input, so their representation requires 
only polynomial space. This characteristic polynomial may be factored into ir- 
reducibles in polynomial time, so the description of p (a) and of all intermediate 
results is polynomially bounded. □ 

It is trivial to check whether a — (3 and whether a belongs to one of N, Z, Q. It 
takes only polynomial time to determine whether a is a root of unity, and if so, 
to calculate its order and phase. 



B Number fields and ideals 

In this section, we define some terminology from algebraic number theory and 
mention textbook results. For more details, see |17|20| . We also define the ideal- 
counting function vp . which is a notion of magnitude of algebraic numbers dis- 
tinct from the usual absolute value. We follow the excellent presentation of |12j . 

An algebraic number field is a field extension K of Q which, considered as a 
Q- vector space, has finite dimension. This dimension is called the degree of the 
number field and is denoted by [K : Q]. The primitive element theorem states 
that for any number field K, there exists an element 9 £ K such that K = 
Q (9). Such a 9 is called a primitive element of K and satisfies ng = [K : Q]. 
The proof of the primitive clement theorem is constructive and shows how to 
obtain a primitive element for K = Q (ai, . . . , ctk) given a%, . . . , afe. There exist 
exactly ng monomorphisms from IK into C, given by 9 — >■ 9i, where 9i are the 



Galois conjugates of 9. If a £ K, then n a \ng. Moreover, if 01, . . . , a ng are the 
mononiorphisms from K into C then u\ (a) , . . . , cr„ e (a) are exactly the Galois 
conjugates of a, each repeated n g /n a times. The norm of a relative to K is 
defined as 

ng 

A/k/q (a) = H<J t (a) = (Afabs (a)) ne/n ° 
»=i 

For a number field K, the set Kz = Az n K forms a ring under the usual 
addition and multiplcation. The ideals of Kz are finitely generated, and form a 
commutative ring under the operations 

IJ = {xy\xe I, y G J} 

I + J = {x + y\x e I,y e J} 

with unit [1] = Kz and zero [0] = {0}. An ideal P is prime if P = AP implies 
A = P or A = [1]. A well-known fact is the fundamental theorem of ideal theory: 
each non-zero ideal may be represented uniquely (up to reordering) as a product 
of prime ideals. 

This theorem gives rise to the following ideal-counting function vp : Kz\ {0} — > 
N. For a fixed prime ideal P, we define vp (a) to be the number of times P ap- 
pears in the factorisation into prime ideals of [a]. That is, 

v P (a) = k if and only if P k \ [a] and P k+1 \ [a] 

We also define vp (0) = oo. The function satisfies the following properties: 

— v P (a/3) = v P (a) + v P (j3) 

— vp (a + (3) > min {vp (a) , vp (/?)} 

— If vp (a) 7^ vp (f3), then v P (a + /?) = min{v P (a) , v P (/3)}. 

For any a € K such that a ^ Kz, we can find (3 s Kz and n € Z C Kz such that 
a = /3/n. We extend wp to IK by defining vp (a) = vp (/?) — vp (n). The first of 
the three properties of vp above guarantees that this value is independent of the 
choice of f3, n, making the extension of vp to K well-defined. Note the extension 
preserves the three properties. 

For an ideal I ^ [0], the quotient ring Kz// is finite. The norm of /, denoted 
N (7), is defined as |K Z /7|. We define also Af ([0]) = oo. Notice that N (I) = 1 if 
and only if I = [1], otherwise Af (I) > 2. Each prime ideal P contains a unique 
prime number p, and Af (P) = pf for some natural number / > 1. In general, 

| A/k/q (a)\ = Af ([<*])> 2^(«) 

Noting that A/" (P) > 2 for any prime ideal P, we have 

v P (a) < log 2 |A/k/ Q (a) | < log 2 \Afabs {a)\ d 

where d = [K : Q]. Thus, if we are given K = Q(ai, . . . for canonically 
represented algebraic numbers Qj and a canonically represented a e K, we can 



observe that d is at most polynomially large in the length of the input and 
\N a bs ( a )\ is a t most exponentially large in the length of the input. Therefore, 
Vp (a) is only polynomially large. 

The following lemma is simple, but occurs frequently in what follows, so we 
state it explicitly here. 

Lemma 2. Let IS, be a number field and a£l with a Kg. Then there exists 
a prime ideal P o/Kz such that vp (a) 7^ 0. 

Proof. There exist fi £ Kz and m £ Z such that a — f3/m. If [fi] — [to], then 
fi and to are associates, so a must be a unit of Kg. Since a ^ Kg, it follows 
that [fi] 7^ [to], so the factorisations of [fi] and [to] into prime ideals must differ. 
Therefore, vp (fi) 7^ vp (to) for some prime ideal P, so vp (a) 7^ 0. □ 



C Baker's theorem and van der Poorten's theorem 

Theorem 1. (Baker 111!) Let a±, ... ,a n , j3%, j3 n be algebraic numbers of de- 
gree at most d. Let Ai > 4 be an upper bound for the height of OLi, and H > 4 be 
an upper bound for the heights of j3\, . . . , /3„ . // 

A = 0i In a,\ + . . . + fin hi a n 7^ 

and 

Q = In Ax . . . In A n 

then 



\xi\A\ > -c\n(HSl)n 



where c = (I6nd) 



2(in„ 



Lemma 3. Let A, b £ C, where \X\ = 1 and X is not a root of unity. Suppose 
(j> (n) is a function from N to C for which there exist a, % G K such that < x < 1 
and |^(n.)| < ax™- There exists an effective bound to suc/i i/iai if 

\ n = 4>{n) + b (12) 

f/ien n < m. Moreover, if A, 6 € A and aj€Q are given as input, then m is at 
most exponential in the length of the input L = ||A|| + + ||a|| + ||x||. 



Proof. The left-hand side of ( 12 ) describes points on the unit circle, whereas the 



right-hand side tends to b. If |6| 7^ 1, then for n large enough, the right-hand 



side of ( 12 1 will always be off the unit circle. This happens when 

ln(||&|-l|/a) 
ln(x) 

The difficult case is when b is on the unit circle. We will use Baker's theorem 
to derive a bound on n. Consider the angle A between A™ and b. This angle can 



be zero for at most one value of n, because A is not a root of unity. Otherwise, 
we have 

A" 

4 = In — = nln(A)-ln& + 2fc„ln(-l) ^0 

where k n is an integer chosen so that A = ir for some r € [0, 27r). Then 2n is 
a height bound for the coefficients in front of the logarithms (because k n < n), 
H = max {H\, H^, 4} is a height bound for the arguments to the logarithms and 
d = max{n>, ?!{,} is a bound on the degrees. Then by Baker's theorem, we have 

In \A\ > - (48d) 600 In 2 H In (2n In 2 H) 
which is equivalent to 

|^|>(2nln 2 i/)- (48!i)60 ° hl2ff 

This is a lower bound on the length of the arc between A™ and b. The length 
of the chord is at least half of the bound: |A™ — b\ > \A\ /2. So in the equation 
\ n — b = 4> (n), the left-hand side is bounded below by an inverse polynomial in n. 
However, the right-hand side shrinks exponentially quickly. For n large enough, 
the right-hand side will forever be smaller in magnitude than the left-hand side. 

We will now quantify the bound on n. Let p\ = (48d) 600 In 2 H and p 2 = 
2 In 2 H. Observe that if A and b are canonically represented algebraic numbers, 



then pi,p2 are polynomials in the size of the input. Then (12 1 cannot hold if 

\ (P2n)- pi > a X n 

which is equivalent to 

- In (2) - In (a) - p\ In (p 2 ) - pi In (n) > n In (x) 

Define p^ = In (2) + In (a) + p\ In (p 2 ) and p^ = max{p3,pi} (also polynomials 
in the size of the input). Then it suffices to have 

p 4 < n 



which is guaranteed by 



In (x) 1 + In (n) 

Vn > 



Mx) 

Observe that — l/ln(x) is at most exponentially large in ||x||- Therefore, the 
bound on n is exponential in the size of the input. □ 

Lemma 4. Suppose A!,A 2 ,a,6,c G C are non-zero, where \Xi\ = |A 2 | = 1 and 
Ai,A 2 are not roots of unity. Let 4>(n) be a function from N to C such that 
< \4> (n)\ < wx n for some w, x € R, X G (0, 1)- Then there exists a computable 
bound m such that if 

a\ , l = b\ 2 l + c + 4>{n) (13) 

then n < m. Moreover, if Ai, A 2 , a, b, c e A and w,x € Q are given, then m is 
at most exponentially large in the length of the input ||Ai|| + ||A 2 || + ||a|| + + 

W + INI + 11x11- ' " " 



Proof. Multiplying the equation by c/ |c| \a\ allows us to assume that \a\ — 1 and 

c e e+. 

Let / (n) — aA™, g (n) — b\ 2 + c. It is clear that / (n) describes points on 
the unit circle 0\ , whilst g (n) describes points on the circle 2 with centre c on 
the real line and radius |6|. 

If these circles do not intersect, then for n large enough, \tp (n) | will be forever 
smaller than the smallest distance between the circles. This happens when 

ln(c-\b\ - 1) -\n(w) 
In (x) 

which is an exponential bound in the size of the input. 

Suppose now the circles intersect in two points, Z\ and z 2 . Let L\ be the 
horizontal line through z\ and Li the horizontal line through z^. Let L\ nOi = 
{£1,2:1}, Li n 02 = {yi,zi}, L 2 n d = {x 2 ,z 2 } and L 2 n 2 = {y 2 ,z 2 }. It is 
trivial that z 2 — ~z\ , X2 = Xi, y 2 = yi- 





We first argue that for n large enough, ( |13| can hold only if for some inter- 
section point Zi, Re(zi) lies between Re(f(n)) and Re(g(n)), or Im(zi) lies 
between 7m (/ (n)) and 7m (5 (n)). This can only be violated in two symmetric 
situations: either / (n) is on the arc ziz 2 of 0\ which lies inside 2 and g (n) is 
on the arc y\y 2 of 2 which lies outside 0\, or / (n) is on the arc X\X 2 of Oi 
which lies outside 2 and <? (n) is on the arc ziz 2 of 2 which lies inside 0\. In 
the first situation, when g (n) is on the arc y\y 2 of 2 outside 0±, we have 

\f(n)~g(n)\ > \g (n)| - 1 > - 1 

This lower bound is positive and independent of n, so equality cannot hold for n 
large enough because <fi (n) tends to zero exponentially quickly. This is the case 
when 

^ ln(| yi |-l)-lnH 



which is exponentially large in the size of the input. The second situation is 
analogous. 

Therefore, we can assume that one of the intersection points z$ separates / (n) 
and g (n) horizontally or vertically. We will show a lower bound on \f (n) — g (n)\ 
which shrinks slower than exponentially. The horizontal and vertical cases are 
completely analogous. We show the working for the horizontal case. Assume that 
Re (zi) lies between Re (/ (n)) and Re (g (n)). Clearly 



1/ (n) - g (n)\ > \Re (g (n) - / (n))| = \Re (z 4 - / (n))\ 
Let a — arg (Ai), 7 = arg (a) and f3 — arg (z,). Then 



\Re (zi - f(n)) \ = |cos(na + 7) - cos(/3)| 



Let u n ,v n be appropriately chosen integers so that 



(3 — na — 7 



j3 + na + 7 



- U n TT G 



' V n TT S 



7T 7T 

2' 2 

7T 7T' 

2' 2 



Then using the inequality 



\Re (g (n) - z t )\ 



ft — na — 7 . /3 + na + 7 



sin ■ 



we have 



|sin (x)\ > — for X G 

7T 



7T 7T 

2' 2 



/3 — na — 7 

sin 

2 

. /3 + na + 7 
sm 



> - 

7T 
1 

> - 
7T 



/3 — na — 7 
2 

/3 + na + 7 



7TM r 



Both of these expressions are sums of logarithms of algebraic numbers, so we 
can give lower bounds for them using Baker's theorem as in Lemma [3] 



\Re( Zi -f(n))\ > (pin) 



-V'l 



for some constants pi,p2 > which are polynomially large in the input. A 
similar lower bound holds for \Re(g(n) — If 8 = arg(A2), f] = arg (6) and 
= arg (z; — c), we have 

|Jfe (3 0) - = H (cos (n<5 + 77) - cos {9)) > (p 3 n)- pi 

where Ps,Pi are positive constants of polynomial size in the input. Hence we 
have 

\f(n)-g(n)\ >2(p 5 n)- p < i 

where p$ = max {pi , j>3 } and pg = max {p2 ,Pi}- Since <fi(n) shrinks exponentially 
quickly, a bound on n follows past which ( 13 1 cannot hold. In the manner of 



Lemma [3] we can show that this bound is exponentially large in the size of the 
input. The vertical case is analogous, except that considering imaginary parts 
gives sines instead of cosines, so we shift all angles by ir/2 and proceed as above. 
If the circles are tangent and neither lies inside the other, then the intersection 
point separates / (n) and g (n) horizontally, so we are done by the above analysis. 

Finally, suppose that the circles are tangent and one lies inside the other: 
|&| + c = 1. The argument of / (n) is 7 + na. By the cosine theorem applied to 
the triangle with vertices / (n) and the centres of the circles, we have 

|/ (n) - c| 2 = c 2 + 1 - 2c cos (7 + na) 

Therefore, the shortest distance from / (n) to a point on O2 is 



h (n) — y/ c 2 + 1 — 2c cos (7 + na) — (1 — c) 

Let A (n) = ^/c 2 + 1 — 2c cos (7 + na) and B = 1 — c. Since A < 1 + c, we have 
A + B < 2, so 

A 2 - B 2 

h (n) = A — B = — ; > c (1 — cos (7 + na)) 

\ i A + B ~ y w " 

Let k n be an integer, so that 

7 + na + fc„27r € [— tt, it) 
A lower bound on this angle follows from Baker's theorem: 

| 7 + na + k n 2n\ > (p 7 n)~ P8 
for some constants pj,p% > which are polynomially large in the input. Then 

cos (7 + na) < cos (^(p7n)~ Pa ^J 

h{n) > c(l-cos((p 7 n) _P8 )) 
From the Taylor expansion of cos (x). it follows easily that 

1 — cos (x) > — x 2 for x < 1 
V ' ~ 24 

Since P7,ps > 1, we have {p7n)~ Ps < 1. Therefore, 

h(n)>c^(p 7 n)- 2p » 

This lower bound on h (n) shrinks inverse-polynomially as n grows. Recall that 
h (n) is the smallest distance from / (n) to O^. It follows that for n large enough, 
\4> {n)\ < h (n) forever, so / (n) = g (n) + <fi (n) cannot hold. In the manner of 
Lemma [3] we can show that the bound is exponentially large in the input. □ 



so 



Theorem 2. (van der Poorten \2V$) Let cti, . . . , a n be algebraic numbers of de- 
gree at most d belonging to a number field K and with heights respectively not 
exceeding A\ , . . . , A n . Let P be a prime ideal of IK containing the rational prime 
p. The following inequalities have no solutions in rational integers b%, . . . , b n with 
absolute values at most B > e 2 : 

oo>v P (ofr . . . ofc - l) > (16 (n + 1) d) 12[n+1) (p d / In (p)) Q (In (B)f 

where fl = In (A\) ... In (A n ). 



D Skolem's problem, depth 2 

In this section, we consider the problem of finding whether a linear recurrent 
sequence u n of depth 2 contains as an element. The characteristic equation 
of the recurrence may have one repeated root 9^0, or two distinct roots #1,2, 
giving either 

u n = (A + Bn) 9 n 

or 

u n = A6 r l + B6™ 

Solving the problem in the former case is trivial. The latter case is an instance 
of the algebraic number power problem: decide whether there exists n g N such 
that 

a n =f3 (14) 

for given a, ft £ A. The algebraic number power problem is decidable |12| . Kan- 
nan and Lipton pQ proved a polynomial bound on n when ft has the form p (a) 
for a given p € Q [x] and a is not a root of unity. We give a brief recapitulation 
of the decidability proof and extract a polynomial bound on n from it. 

Lemma 5. Suppose a, ft £ A. If a is not a root of unity, then there exists 



a computable bound m such that if (14) holds, then n < m. Moreover, m is 



polynomial in the length of the input \\a\\ + \\ft\\- 

Proof. Let K — Q (a, ft). If a is not an algebraic integer, then by Lemma[2]there 
exists a prime ideal P in the ring such that vp (a) 7^ 0. Then if a n = ft, we 
have 

vp (a n ) — nvp (a) — vp (ft) 
If vp (a) and vp (ft) have different signs, then we are done. Otherwise, 

n = < \v P G8)| < log 2 |A/k/q GS)| < log 2 \M abs (ft)\ d 

where d = [Q (a, ft) : Q] is at most polynomially large in ||a|| + \\ft\\. It follows 
that the bound on n is polynomially large in the length of the input. 



Suppose a is an algebraic integer. It is not a root of unity, so by Kronecker's 
theorem |22J, a has a Galois conjugate a (a) with magnitude strictly greater than 
1. In fact, a significant strengthening of Kronecker's theorem, due to Blanksby 
and Montgomery [53], guarantees the existence of a conjugate a (a) such that 

\a(a)\>l + 



30n^ In (6n a ) 
which implies 

j — | — 7— r; < 60n^ In (6n Q ) 
ln|cr(o;)| 



Then if a™ = 0, we have 

ln|fT(/?)l <m|a(/3)|60n>(6n c 



i»k(«)l 

Observe that if we are given a canonical description of (3, both ln|cr(/3)| is at 
most polynomially large in ||/3||, whereas 60n^ln(6n Q ) is at most polynomially 
large in ||a||. It follows that the bound on n is polynomial in the length of the 
input. □ 

The condition that a is not a root of unity is obviously necessary in Lemma 
[5] because if j3 is also a root of unity, a n — (3 could hold for all n in a linear set. 
It is easy to exhibit linear recurrences of depth 2 with infinitely many zeroes (for 
example, u\ = 0, 112 = 1, u n +2 = u n ), but by the Lemma, they all have have 
two distinct roots whose ratio is a root of unity. 



E Skolem's problem, depth 3 

In this section we will focus on Skolem's problem for linear recurrent sequences 
of depth 3. The characteristic equation of such a sequence may have either three 
distinct roots a, j3, 7, or one repeated root a and one simple root /?. If the three 
roots are distinct, we are concerned with solving for n € N equations of the form 

Aa n + Bf3 n + C 7 ™ = (15) 

where A, B, C, a, f3, 7 S A are given and non-zero (if any of them is 0, then the 



sequence satisfies a recurrence relation of smaller depth). Then ( 15 1 is equivalent 
to 

'/3\™ C n\ n A 



B KaJ B ^ 

We will consider only non-singular sequences: the ratios of the roots a,/?, 7 are 
not roots of unity. Let also |a| > \/3\ > | 7 |. In Lemmas [6j [7] [8] below, the length 
of the input is || A|| + ||B|| + ||C|| + |M| + \\/3\\ + || 7 ||. 

Lemma 6. // |a| > \0\, then there exists an effective bound m such that if 



equation (15) holds, then n < m. Moreover, m is at most exponential in the 



length of input. 



Proof. Trivial, because Aa" is dominant. If 

'1n\A/2B\ ln\A/2C\\ 



n > max 



In 



\p/a\ ' ln| 7 /a| J 



then 



□ 



A\a A\aJ 



< 



A la 



c ny 

A \a) 



1 1 

<2 + 2 =1 



3\ > |7|, then there exists an effective bound m such that 
if equation (15) holds, then n < m. Moreover, m is at most exponential in the 



Lemma 7. If \a\ 



length of the input. 

Proof. This is a direct application of Lemma [3] to equation ( 16 1 . □ 



Lemma 8. // \a\ 



3\ = I7I, then there exist at most two values of n such that 
equation (15) holds. Moreover, they are at most exponential in the length of the 



input and are computable in polynomial time. 



Proof. The left-hand side of (16) as a function of n describes points on the unit 



circle in the complex plane, whereas the right-hand side describes points on a 
circle centred at — A/B with radius \C/B\. Note these circles do not coincide, 
because A ^ 0. We can obtain their equations and compute their intersection 



point(s). If they do not intersect, then equation (15 1 can never hold. Otherwise, 
the equation can only hold if the two sides are simultaneously equal to the same 
intersection point. For each of the (at most two) intersection points 8, let 



Si 



S-2 



c ny 

B W 



A 
B 



Observe that \Si\ < 1, because (3/a and 7/0 are not roots of unity. We compute 
Si and S2 from the bound in Lemma [5] and check whether Si is non-empty. 
□ 

Next, we consider recurrent sequences of depth 3 with one repeated and one 
simple root. We are given non-zero A, B, C,a,(3 £ A, and the length of the input 
is || A|| + ||B|| + ||C|| + ||a|| + \\/3\\. We wish to solve for n 



(A + Bn) a n + Cfi n = 



(17) 



Lemma 9. There exists an effective bound m such that if holds, then n < m. 
Moreover, m is at most exponential in the length of the input. 



Proof. If | a | > |/3 1, then for 



n > 



\A\ + \C\ 
\B\ 



we have 



so 



|C| < \B\n-\A\ < \A + Bn\ 



\Cp n \ < \(A + Bn)a n \ 



and 17 cannot hold. Now suppose \a\ > \j3\ and rewrite (17 1 as 



A + Bn 
C 



A bound on n follows from 



which is implied by 



n < 



d(n + 1) < 



where d = max {\A/C\ ,\B/C\}. Taking logarithms, we see that it suffices to 
have 

n f 
l + ln(n + l) > m(/3/o) 

where / = max {In (d) , 1}. Noting that 1 + In (n + 1) < lyjn for all n > 1, we 
see that it suffices to have 

n>4/ 2 /ln 2 {fi/a) 
This is an exponential bound in the length of the input. □ 



F Skolem's problem, depth 4 

In this section we will give lemmas which form a decidability proof for Skolem's 
problem for depth 4. Assume algebraic numbers A,B,C,D and a, /3,7,<5 are 
given and the input has length ||yl|| + ||S|| + ||C*|| + \\D\\ + \\a\\ + \\/3\\ + || 7 || + ||<J||. 
We wish to solve for n the following equations: 



Aa n 


+ Bf3 n + C* 7 ™ + DS n = 


= 


(18) 


{A + 


Bn) a n + CP' n + D~/ n 


= 


(19) 


{A + 


Bn) a n + {C + Dn) /3 n 


= 


(20) 


(A 4 


Bn + Cn 2 ) a n + D/3 n 


= 


(21) 


(A 4 


Bn + Cn 2 + Dn 3 ) a n 


= 


(22) 



As before, we assume that the ratios of a, f3, 7, <5 are not roots of unity. Solving 
(22) is trivial. We can rearrange (21 1 as 



(A + Bn + Cn 2 



-D 



(23) 



The left-hand side tends to or 00 in magnitude, depending on whether |a| < |/3| 
or not. In both cases, a computable bound on n follows. The remaining equations 
(18l)pjlp0l are more involved. 



Lemma 10. There exists an effective bound m such that if equation (20) holds, 
then n < m. Moreover, this bound is at most exponential in the size of the input. 

(C 



Proof. Rearrange ( 20 1 as 



A" = -- 



Dn) 



(24) 



(A + Bn) 

where A = a/ ft- The right-hand side of p4| tends to —D/B as n tends to infinity. 

First, if A is an algebraic integer, then by Blanksby and Montgomery's the- 
orem [23 , it has a Galois conjugate a (A) such that 

k(A)| > 1+ 1 



30n 2 l n (6n A ) 



We apply the monomorphism to both sides of (|24j) , so that the right-hand side 
tends to a (D/B) whereas the left increases in magnitude exponentially quickly, 
giving a bound on n. 

Second, suppose A is not an algebraic integer. Then by Lemma[2]there exists 
a prime ideal P in the ring of integers of K = Q (a, f3, A, B, C, D) such that 
vp (A) 7^ 0. Without loss of generality, we can assume vp (A) > (if vp (A) < 0, 
swap a with f3, A with C, and B with D). Applying vp to (24 1 gives 



vp (A™) = nvp (A) =vp[ — 



C + Dn 
A + Bn 



< ln 



■A/k 



K/Q 



( C+Dn \ 
\ A+Bn J 



C + Dn 
A + Bn 

< p\, where e\ 



We will now show that for n > e\ we have ln 
and pi are an exponential and a polynomial constant in the input, respectively. 
This suffices because nvp (A) > n grows with n. 
Consider 

, C + Dn 
fin) 

By the triangle inequality, we have 
\D\n-\C\ 



A + Bn 



h(n) 



<f(n) < 



\D\n- 
\B\n 



\C\ 
\A\ 



= 9{n) 



\B\n + \A\ 

Assume that n > \A/B\. It is easy to see that h and g tend to \D/B\, h from 
below and g from above. For 

2\BC\ +2\AD\ + \AB\ 



n > 



w\ 



we have 



Similarly, for 



we have 



n > 



g(n) < 
2\BC\+2\AD\ - \AB\ 



h(n) > 



\B\ 

D 
B 



Thus, for n at most exponentially large, we have / (n) < \D/B\ + 1. This bound 
is at most exponential in the size of the input. 

Applying the monomorphisms a from K into C to 



gives terms of the form 



C + Dn 
A + Bn 



a{C)+a (D) n 
' a (A) + a (B) n 



We can obtain an upper bound, exponential in the input and independent of n, 
on the magnitude of each such term. There are exactly [K : Q] such terms, one 
for each monomorphism a of K into C. This gives a polynomially large upper 



bound on In 



\r , ( C+Dn \ 



which holds for exponentially large n. □ 



Lemma 11. There exists an effective bound m such that if equation {19) holds, 
then n < m. This bound is at most exponential in the size of the input. 

Proof. First suppose |a| > \/3\ , |7|. Then the term (A + Bn)a n is dominant. 
More precisely, rewrite (19) as 



A + Bn = -C 



a 



-D 



and observe that if 



n > 



\A\ + \C\ + \D\ 
\B\ 



then 



\A + Bn\ >\B\n- \A\ > \C\ + \D\ > 



-C 



a 



-D 



so (19) cannot hold. 



Second, suppose that \f3\ > \a\ , Then the term C/3™ is dominant. More 



precisely, rewrite ( 19 1 as 

(A + Bn) 



D 



= -C 



We will require that 



and 



D 



< 



\C\ 



(A + Bn) ( - 



< 



\C\ 



The former inequality holds for n > In \C/2D\ / In | T / /? | , which is at most expo- 
nentially large in the input. The latter inequality is implied by 



(n + 1) 



< 



\_C\_ 
2M 



where M = max{|A| , Now let r = [-In (2) /\n{a//3)~\, so that 

1 



< 



and consider only n of the form n = kr for k G Z + . If 

ln|C/4Mr| 



k > 



In (7/8) 



and k > 5, we have 



k< ( - ) (fe + l)< ( - 



< 



4Mr 



so 



(n+1) < 



2n < 



2M 



It is clear that r is at most exponentially large in the size of the input, whereas 
the bound on k is polynomial. Therefore, the bound on n is exponential. 
Finally, suppose |/3| = I7I > \a\. Rewrite (19 1 as 



1, 



D A + Bn ( a\ n 
C ' 



c 



Then an exponential bound on n follows from Lemma [3] because the right-hand 
side is a constant plus an exponentially decaying term, whereas the left-hand 
side is on unit circle. □ 



Lemma 12. //a,/3,7, S do not all have the same magnitude, then there exists 
an effective bound m such that if equation (18) holds, then n < m. Moreover, m 
is at most exponentially large in the input size. 



Proof. Let \a\ > \/3\ > | 7 | > \S\. First, if \a\ > then Aa n is the dominant 



term in (18 1. Rewrite the equation as 



D 
A \a 



c ny 

A \aJ 



A [a 



and observe that if 



n > max 



]n\3B/A\ ln|3C/A| \n\3D/A\ \ 
ln|a//3| ' ln|a/ 7 | ' ln|a/J| J 



then 



C 
A 



(IT 



D 
A 



Second, if |a| = > I7I, then rewrite (18) as 



.4 
B 



c n\ n 

B \cu 



< 



D 
B 



(25) 



The left-hand side of (25 1 is on the unit circle, whereas the right is a constant plus 



exponentially decaying terms. An exponential bound on n follows from Lemma 

El 

Finally, if |a| — \f3\ = \ j\ > \S\, then an exponential bound on n follows from 
Lemma |4] applied to equation ( p5| . □ 



Thus, the only outstanding problem is to solve (18 1 when |or| = \/3\ — \j\ — 
\S\. This is difficult for general a,/?,7, 8, so we will restrict ourselves to two 
sufficient special cases: when at least two of them are real, and when they are 
two pairs of complex conjugates. We will also assume that they are algebraic 
integers. This is sufficient for our application to the orbit problem, because any 
orbit instance 3n.A n x € V may be reduced in polynomial time to an instance 
where A is an integer matrix, so that its eigenvalues are algebraic integers. This 
trick is not readily available when solving a Skolem instance. 

Lemma 13. Let a, /3,7,<5 be algebraic integers with \a\ = = I7I = \S\. If 
a,/3 € K or (a,/3) and (7,^) are pairwise complex conjugates, there exists an 
effective bound m such that if equation (18) holds, then n < m. Moreover, m is 
exponential in the length of the input. 



Proof. Let K = Q (a, /3, 7, 5, A, B, C, D). First suppose that a, f3 e 
then we have a Skolem instance of depth 3 

3n. (A + B)a n + C^ n + Dj n = 



If a = f3, 



An exponential bound on n follows from Lemma [8] If a = — j3, then we consider 
even n and odd n separately, obtaining two Skolem instances of depth 3 of the 
same type. 

Now suppose j3 = a and 7 = 5. If a/ f3 is an algebraic integer, then since it 
is not a root of unity, there exists a monomorphism a from K to C such that 



\a (a) | 7^ | cr (/3)|. Applying it to ( 18 1 leads to a Skolem instance of depth 4 with 
roots whose magnitudes are not all the same. A bound on n follows from Lemma 

m 

Suppose then that a/ (3 is not an algebraic integer. By the reasoning of Lemma 
[2] there exists a prime ideal P in Kz such that vp (a) 7^ vp (/3) and at least one 
of vp (a) and vp (/?) is strictly positive. Assume without loss of generality that 

v P (a) > v P 03) > 

Since a/3 — jS — \a\ 2 , we have 

v P (a) + v P {fi) = v P (7) + v P (6) 

Therefore, at most two of the roots are smallest under the valuation vp. 

If one root, say /3, is strictly smaller under vp than the rest, then rewrite 



(18j) as 

Aa n + BP 11 = -Cj n - DS n (26) 
Since vp (j3) < vp (a), for n > vp (A/B) jvp (/?/«) we have 
v P {Aa 11 + Bfi n ) = v P (B) + nv P (f3) 

whereas 

v P (-C7™ - DS n ) > v P (C) + nv P (7) 



Therefore, for n > vp (B/C) jvp (7//?), we have that the left-hand side of (26 1 



is strictly smaller under vp than the right-hand side, so (18 1 cannot hold. This 
bound on n is polynomial in the input size. 

Now suppose that there are two roots with strictly smallest valuation: 

< vp (j3) = vp (7) < vp (a) = vp (S) 



Then rewrite ( 18 1 as 



Bp n ( ( ] ( 1 ) - 1 ) = Aa n + D5 n (27) 



B ) \(3 

Since 7//? is not a root of unity, the term (-C/B) (7//?)™ — 1 can be zero for 
at most one value of n. This value is at most polynomially large in the input 
size (by Lemma [5|. For all other n, we may apply van der Poorten's theorem 
to this term. Let p be the unique prime rational integer in the ideal P, and let 
d = [K : Q]. Let H be an upper bound for the heights of —C/B and 7//?. Then 
by Theorem [2j we have 

-(B) (i)"- i ) s(4M,36 4"" ff)2<i,i " )2 

It is classical N (P) — p* for some positive integer /, so A/"(P) > p. Moreover, 
since a is an algebraic integer, all prime ideals Pi, . . . , P s in the factorisation of 
[a] appear with positive exponents fci, . . . , k s 4 - 

[a]=Pt...P*° 



Since TV (Pi) > 2 for all Pi, we have 

I/Vk/q (a)\ =Af([a\) >N(P)>p 



Therefore, p is at most exponentially large in the input size. Then we can write 
the inequality from van der Poorten's theorem as 



v P 



- lj < Ex (Inn) 

where E\ is exponentially large in the input size and independent of n. Now we 
apply vp to both sides of equation (27 1: 

v P (LHS) < v P (B) + nv P {fi) + E l (Inn) 2 



and 



v P (RHS) > v P (A) + nv P {a) 



Equation (18 I cannot hold if 

v P (B) + nv P (/?) + Ei (Inn) 2 < v P (A) + nv P (a) 
which is implied by 

v P (B/A) + E 1 (Inn) 2 <n 
since vp (a) > vp (/?). Let E% = max{?ip (B/A) ,Ei}, then this is implied by 



Since 



E 2 ((Inn) + lj < n 

5-y/n 



(lnn) z + 1 < 



for all n > 1, it suffices to have 



n > I -E 2 



This bound on n is exponential in the size of the input. □ 



References 

1. Kannan, R., Lipton, R.: Polynomial-time algorithm for the orbit problem. Journal 
of the ACM 33(4) (1986) 808-821 

2. Arvind, V., Vfjayaraghavan, T.: The orbit problem is in the GapL hierarchy. J. 
Comb. Optim. 21(1) (2011) 124-137 

3. Tarasov, S.P., Vyalyi, M.N.: Orbits of linear maps and regular languages. CoRR 
abs/1011.1842 (2010) 



4. Skolem, T.: Ein Verfahren zur Behandlung gewisser exponentialer Gleichungen 
und diophantischer Gleichungen. Skand. Mat. Kongr. 8 (1934) 163-188 

5. Mahler, K.: Eine arithmetische Eigenschaft der Taylor-Koeffizientcn rationaler 
Funktionen. Proc. Akad. Wet. Amsterdam 38 (1935) 51-60 

6. Lech, C: A note on recurring series. Arkiv for Matematik 2 (1953) 417-421 

7. Hansel, G.: Une demonstration simple du theoreme de Skolem-Mahler-Lech. The- 
oretical Computer Science 43 (1986) 91 - 98 

8. Halava, V.: Decidable and undecidable problems in matrix theory. TUCS Technical 
Report 127 (1997) 

9. Vereshchagin, N.: Occurrence of zero in a linear recursive sequence. Mathematical 
Notes 38 (1985) 609-615 

10. Mignotte, M., Shorey, T., Tijdeman, R.: The distance between terms of an algebraic 
recurrence sequence. Jour. Reine Angew. Math. 349 (1984) 63 - 76 

11. Baker, A.: Transcendental number theory. Cambridge University Press (1975) 

12. Halava, V., Harju, T., Hirvensalo, M., Karhumaki, J.: Skolem's problem - on the 
border between decidability and undecidability. TUCS Technical Report (683) 
(2005) 

13. Blondel, V., Portier, N.: The presence of a zero in an integer linear recurrent 
sequence is NP-hard to decide (2002) 

14. Schonhage, A.: On the power of random access machines. In Maurer, H., ed.: 
Automata, Languages and Programming. Volume 71 of Lecture Notes in Computer 
Science. Springer Berlin / Heidelberg (1979) 520-529 

15. Mignotte, M.: Some useful bounds. Computer Algebra (1982) 259-263 

16. Pan, V.: Optimal and nearly optimal algorithms for approximating polynomial 
zeros. Computers & Mathematics with Applications 31(12) (1996) 97 - 138 

17. Cohen, H.: A Course in Computational Algebraic Number Theory. Springer (1993) 

18. Lenstra, A.K., Lenstra, H.W., Lovasz, L.: Factoring polynomials with rational 
coefficients. Mathematische Annalen 261 (1982) 515-534 

19. Zippel, R.: Zero testing of algebraic functions. Information Processing Letters 
61(2) (1997) 63 - 67 

20. Stewart, I., Tall, D.: Algebraic Number Theory and Fermat's Last Theorem. 3rd 
edn. A. K. Peters (2002) 

21. van der Poorten, A.J.: Linear forms in logarithms in the p-adic case. Transcendence 
Theory: Advances and Applications (1977) 29-57 

22. Kronecker, L.: Zwei Satze iiber Gleichungen mit ganzzahligen Koeffizienten. J. 
Reine Angew. Math. 53 (1875) 173-175 

23. Blanksby, P., Montgomery, H.: Algebraic integers near the unit circle. Acta Arith. 
(1971) 355-369 

24. Shidlovskii, A.: Transcendental Numbers. New York: de Gruyter studies in math- 
ematics (1989) 

25. Kozen, D.: The Design and Analysis of Algorithms. Springer- Verlag, New York 
(1991) 

26. Chonev, V.: The orbit problem in zero and one dimensions. Master's thesis, Oxford 
University Department of Computer Science (2011) 



